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^ ; ABSTRACT 

<~i , In this paper we report the resuhs of high-resolution circular spectropolarimetric monitoring 

Qh' of the Herbig Ae star V380 Ori, in which we discovered a magnetic field in 2005. A careful 

study of the intensity spectrum reveals the presence of a cool spectroscopic companion. By 
Jl^ . modelling the binary spectrum we infer the effective temperature of both stars: 10500 + 500 K 

^ ' for the primary, and 5500 + 500 K for the secondary, and we argue that the high metallicity 

^ \ ([M///] = 0.5), required to fit the lines may imply that the primary is a chemically peculiar 

star We observe that the radial velocity of the secondary's lines varies with time, while that 
of the the primary does not. By fitting these variations we derive the orbital parameters of the 
system. We find an orbital period of 104 + 5 d, and a mass ratio (Mp/Ms) larger than 2.9. 
■ The intensity spectrum is heavily contaminated with strong, broad and variable emission. A 

simple analysis of these lines reveals that a disk might surround the binary, and that a wind 
occurs in the environment of the system. Finally, we performed a magnetic analysis using the 
• Least-Squares Deconvolved (LSD) profiles of the Stokes V spectra of both stars, and adopting 

the oblique rotator model. From rotational modulation of the primary's Stokes V signatures, 
^•—•^ i we infer its rotation period P = 4.31276 + 0.00042 d, and find that it hosts a centred dipole 

' magnetic field of polar strength 2.12 + 0.15 kG, with a magnetic obliquity /3 - 66 + 5°, and a 

I rotation axis inclination / = 32 + 5°. However, no magnetic field is detected in the secondary, 

and if it hosts a dipolar magnetic field, its strength must be below about 500 G, to be consistent 
with our observations. 
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1 INTRODUCTION 

Today, one of the greatest challenges in stellar physics is to under- 
stand the origin of the magnetic fields observed in stars across the 
H-R diagram. It is now well established that the magnetic fields of 
the sun and other cool, low-mass stars are generated by a convective 
dynamo occurring in their external envelope. Their magnetic fields 
are of complex structure, highly variable on short times-scales, and 
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(CFHT) which is operated by the National Research Council of Canada, 
the Institut National des Sciences de TUnivers of the Centre National de la 
Recherche Scientifique of France, and the University of Hawaii 
t Based on observations obtained at the Bernard Lyot Telescope (TBL, Pic 
du Midi, France) of the Midi-Pyrenees Observatory, which is operated by 
the Institut National des Sciences de TUnivers of the Centre National de la 
Recherche Scientifique of France, 
t E-mail: evelyne.alecian@rmc.ca 



are strongly correlated with fundamental stellar parameters such as 
the rotation period, mass , and age, consistent with the dynamo the- 
ory ( iBaliunas et all 1 9951 : [Saailll 9961 : iDonati et alj2003ll2008l) . 

Around 5% of the main-sequence, intermediate-mass A and B 
stars host organised magnetic fields, with strengths b etween 300 G 
and 3 kG, that are stable over tens of years (e.g. lAuriere et al.l 
l2007h . These characteristics show no correlation with the mass, age 
or rotation period of the stars. These qualitative differences suggest 
a different origin of their magnetic fields. Some authors have pro- 
posed the generation of magnetic fields inside the convective core 
of the intermediate mass stars, with a similar dynamo process as in 
the low-mass stars, that would diffuse through their radiative enve- 
lope to reach the surface. However, this mechanism is not able to 
reproduce the observed field characteristics , nor can it explain the 
long-term stability of these magnetic fields ( iMosslllOOlh . 

Currently, the favoured scenario for the origin of magnetic 
fields in A and B type stars is the so-called "fossil field" hypoth- 
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esis. This model proposes that the magnetic fields observed in the 
intermediate mass stars are relics, either of the Galactic magnetic 
field present in the molecular clouds from which the stars formed, 
or generated by a dynamo during the first stages of stellar forma- 
tion. Once a star is bom, the relic magnetic field would necessarily 
survive the various stellar evolutionary stages from the pre-main 
sequence (PMS) to the post-main sequence phases, without signif- 
icant regeneration. 

Magnetic fields are observed among molecular clo uds with 
magnetic strengths between 1 and 100 fiG ( lHeiles|[l997l) and also 
among the main sequence (MS) A/B stars. Until recently we were 
missing information about the magnetic properties of the PMS in- 
termediate mass stars. According to the fossil theory 1-10% of PMS 
intermediate mass stars should be magnetic with similar magnetic 
topology than in the MS A/B stars, and with magnetic strengths 
consistent with the conservation of the magnetic flux from the PMS 
to the MS stages. 

The Herbig Ae/Be stars have been defined bv lHerbij ( Il960f) as 
emission line stars of spectral type A or B, associated with nebulae, 
and situated in an obscured region. These characteristics strongly 
suggest that the Herbig Ae/Be (HAeBe) stars are very young, still 
completing their PMS phase. We therefore consider that the HAeBe 
stars are the PMS progenitors of the MS A/B stars. Many au- 
thors have predicted the presence of magnetic fields in the se stars 
dCatala et al J 1 19991 : iHubrig et all l2004l : IWade et al.l l2007h . Some 
of them try to detect them , without much success, except in HD 
104237 toonatietal.lll997h . 

In order to provide new observational constraints on the fossil 
field hypothesis and explore the magnetic properties of the PMS 
intermediate mass stars, we performed a survey of ~ 70 HAeBe 
stars with the new generation of high-resolution spectropolarime- 
ters: ESPaDOnS, installed at the Canada-France-Hawaii Telescope 
(CFHT) in Hawaii, and Narval, installed at the Bernard Lyot Tele- 
scope (TBL) at Pic du Midi in France (Alecian et al. in prep.). 
Among this sample, we have disco vere d 4 ne w magneti c HAeBe 
stars, including V3 80 Ori jWade etaf] boOSi : ICatalaet al.l l2007l : 
lAlecian et al. l l2008l : lFolsometalJ I l2008h . Our survey has thereby 
brought the first strong arguments in favour of the fossil theory. In 
particular, we have concluded that ~ 6% of the Herbig Ae/Be stars 
observed in our survey are magnetic. 

To go further, we have characterised the magnetic fields of 3 of 
these stars (HD 200775, Alecian et al. 2008; HD 190073, Catala et 
al. 2007; HD 72106, Folsom et al. 2008) in order to compare their 
magnetic topologies and strengths with those of the MS A/B stars. 
We have found that the magnetic Herbig Ae/Be stars host large- 
scale organised magnetic fields, with important dipole components 
with polar strengths between 300 G and 1 .4 kG, in agreement with 
those of the main sequence stars and with the fossil model. 

V380 Ori, of spectral type B9, has been classified as Her- 
big Ae/Be star by Herbig (1960). This star is associated with 
many Herbig-Haro objects , with the Hii rep on NGC 1999 in the 
Orion OBlc association teorcoran&Rav 1995^, and its spec- 
trum contains str ong emission lines ( Finkenzeller & Mundt 1984|; 
iFinkenzeller & Jankovicslll984h . [Hillenbrand et alj ( 119921) classify 
V380 Ori as a Group I object whose spectral energy distribution 
(SED) in the infra-red (IR) is similar to those of the low-mass PMS 
T Tauri stars. The SED is well reproduced with a model of a flat, op- 
tically thick accretion disk. These results firmly establish the PMS 
status of V380 Ori. 

V38 Ori was detected as a rath er bright X-ray source by 
ROSAT ( IZinnecker&Preibisch|ll994h. Subsequent X -rav grat- 
ing spectroscopy with Chandra ( Hamaguchi et aljr2005h implied a 



single-temperature plasma with a relatively high temperature, lead- 
ing these authors to speculate that magnetic activity was occurring 
in the close enviror mient of the star. However, as pointed out by 
IStelzer et alj ( l2005t) . due to the small binary separation of V380 
Ori (see Sect. 3.1) Chandra was not able to resolve the system, and 
the origin of the (apparently single) X-ray source remains obscure. 

The magnetic field of V380 Ori was detected dur ing the first 
run of our Herbig survey, using ESPaDOnS at CFHT ( Wade et alj 
2005). This paper reports subsequent observations of V380 Ori ac- 
quired between 2005 and 2007, as well as the analysis that we 
have performed in order to characterise its magnetic field. In Sec- 
tion [2] we describe the observations and data reduction. Section [3] 
describes the intensity spectrum and reports the discovery of a spec- 
troscopic companion. In Section |4] we determine the fundamental 
properties of both components, and in Section [5] we describe the 
magnetic field analysis performed on both the primary and sec- 
ondary stars. Finally a discussion and a conclusion are given in 
Section|6] 



2 OBSERVATIONS AND DATA REDUCTION 

Our data were obtained using the high resolution spectropolarime- 
ters ESPaDOnS, installed on the 3.6 m Canada-France-Hawaii 
Telescope, and Narval, installed on the 1.9 m Bernard Lyot Tele- 
scope (Donati et al., in preparation) during many observing runs 
spread over nearly 3 years from Feb. 2005 to Nov. 2007. We ob- 
tained 24 circular polarisation (Stokes V) spectra covering the en- 
tire optical spectrum from 3690 to 10480 A, with a resolving power 
/i/A/l = 65,000, and with signal to noise (S/N) ratio between 70 
and 290. Table[T]presents the log of the observations. 

We used both instruments in polarimetric mode. Each Stokes 
V observation was divided into 4 sub-exposures of equal time in 
order to compute the optimal extraction of the polarisation spectra 
(Donati et al. 1997, Donati et al., in preparation). We recorded only 
circular polarisation, as the Zeeman signature expected in linear 
polarisation is about one order of magnitude weaker than circu- 
lar polarisation. The data were reduced using the "Libre ESpRIT" 
package especially developed for ESPaDOnS and Narval, and in- 
stalled at the CFHT and at the TBL (Donati et al. 1997, Donati 
et al., in preparation). After reduction, we obtained the intensity 
Stokes / and the circular polarisation Stokes V spectra of the star 
observed, both normalised to the continuum of V380 Ori. A null 
spectrum (N) is also computed in order to diagnose spurious polar- 
isation signatures, and to provide an additional verific ation that the 
signatures in the Stokes V spectrum are indeed real dDonati et al.l 



3 PROPERTIES OF THE INTENSITY SPECTRUM 
3.1 V380 Ori: an SB2 system 

Many authors have mentioned the pres ence of a companio n around 
the variable primary star of V380 Ori. Ijonckheerd ( Il917h was the 
first to observe a companion of visual magnitude equal to 13 mag 
at a distance of 3" and a position angle of 220°. However these 
observations were not confirmed by Herbig (1960), who proposed 
that Jonckhere was observing bright nebulosities in the vicinity of 

the central star. 

Later, iLeinert et al.l ( Il994ll997h detected an IR companion 
using speckle-interferometry, a detection that was confirmed by 
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Table 1. Log of the observations. Columns 1 and 2 give tlie UT date and the Heliocentric Julian Date of the observations. Column 3 gives the total exposure 
time. Column 4 gives the peak signal to noise ratio (per ccd pixel, integrated along the direction perpendicular to the dispersion) in the spectra, at the wavelength 
indicated in column 5. Columns 6 and 7 give the signal to noise ratio in the reconstructed LSD Stokes V profiles of the primary (P) and the secondary (S). 
Column 8 gives the longitudinal magnetic field of the primary, column 9 gives the rotation phase of the primary as derived in Section 5.2, and columns 10 and 
1 1 give the radial velocities of both components of the system. Column 1 2 gives the instrument used. 
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ISmith et alj (l2005h and iBaines etaP j200d ) using the same tech- 
nique. The angular separation between the primary an d the IR com- 
panion is 0.154", ie 62 AU at a distance of 400 pc jLeinert et al.l 
1 1 997h . Ib ailevi ( 1 1 998h obtained spectroscopic observations of a sam- 
ple of PMS stars, including V380 Ori, 2 years after the obser- 
vations of Leinert (1994,1997), and detected a companion using 
the spectro-astrometry technique. They measured a position angle 
of 207°, which is very close to the measurement of Leinert et al. 
(1997), and they suggested that the fainter c omponent of the sys- 
tem h as the greater equivalent width of Hor. iBouvier & Corporonl 
( I2OOI ') analysed the system using high-angular resolution spectro- 
imaging. They detected the companion in December 1993 in only 
the 3 reddest photometric bands J, H and K. They derived a separa- 
tion between the primary and the secondary of 0. 14", and a position 
angle of 205°. All of these results are generally consistent with the 
presence of an IR-bright but optically-faint companion star sepa- 
rate d from the B9 pri mary by ~ 0.15". 

ICorcorarJ jl994h detected the Li i 6707 A line in the spec- 
trum of V380 Ori, confirmed later by Corporon & Lagrange' ('l999|) 
and Blondel & Tjin A Djie (2006). Corporon & Lagrange ( 199%) 
attributed the Li line to a companion of spectral type later than F5, 
but they did not detect any radial velocity variations in the spectral 
hnes of the primary. 

In our spectra we confirm the detection of the Li 1 6707 A line, 
and we detect the presence of other lines belonging to a cool com- 
panion (Fig. [TJ. We observe that the radial velocity of the lines of 
the cool companion (hereinafter the "secondary"), including the Li i 
line, vary with a time scale of a month, while the radial velocity of 
lines of the B9 HAeBe star (the "primary") do not vary detectably 



with time (see Fig.[TJ. The study of the radial velocities presented 
in Sect. 4 leads to an estimate of the upper limit of the separation 
between both stars: a sin / < 0.33 AU, which is around 200 times 
smaller than the separati on observed b etwee n the IR companion 
and the primary star by iLeinert et alj jl997h . For these reasons, 
we conclude that the spectroscopic companion cannot be the IR- 
bright companion observed using speckle-interferometry, and that 
V380 Ori is therefore a triple system. In the following the primary 
(V380 Ori A - the B9 HAeBe star) and the secondary (V380 Ori 
B - the cooler companion) stars refer to the spectroscopic compo- 
nents, the primary being the hottest and most luminous. The ter- 
tiary refers to the IR-bright companion. The large luminosity ra- 
tio of the spectroscopic components, their small angular separation 
(a sin; < 0.83 mas, assuming a distance of 400 pc), and the rela- 
tively faint visual magnitude {V = 10.2) of the system, explain why 
the secondary has never been directly detected in the past. 



3.2 Fitting of the intensity spectrum 

In order to determine the effective temperatures of both stars and 
their luminosity ratio, we have fitted the intensity spectrum of V380 
Ori with a synthetic spectrum calculated as follows. We used the 
code BINMAGl developed by Oleg Kochukhov to calculate the 
composite spectrum of the binary star. This code takes as input two 
synthetic spectra of different effective temperatures and gravities, 
each corresponding to one of the two components. The code con- 
volves the synthetic spectra with instrumental, turbulent and ro- 
tational broadening profiles, and combines them according to the 
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Figure 1. Selected regions of the spectrum of V380 Ori observed on Jan- 
uary 9th (black line) and February 1 1th (grey line) 2006 around 6425 A (up) 
and around 6400 A (bottom). The positions and ions of the lines are indi- 
cated with a vertical black line for the January spectrum and with a vertical 
grey line for the February spectrum. The labels P and S beside the ion iden- 
tification indicate the primary and secondary, respectively. Note the radial 
velocity variation of the lines of the secondary component, while no shift is 
observed in the radial velocity of the primary's lines. 



radii ratio of the components specified by the user, and the flux ra- 
tio at the considered wavelength given by the atmosphere models, 
to produce the spectrum of the binary star. The individual synthetic 
spectra were calculated in the local thermodyn amic equilibrium 
(LTE) approximation, using the code SYNTH o f lPiskunovl l fT993) . 
SYNTH requires, as inpu t, at mosphere models, obtained using the 
ATLAS 9 code ( lKuruczlll993l) . and a list of spectral line data ob- 
tained from the VALD databas43 (Vienna Atomic Line Database). 

We compared the observed spectrum of V380 Ori to synthetic 
binary spectra using effective temperatures ranging from 9500 to 
11500 K for the primary and from 4500 to 6500 K for the sec- 
ondary, and varying the ratio of radii from 1.0 to 1.8. In order to 
determine the temperature of the primary, we first focussed on re- 
gions of the spectrum containing no significant circumstellar emis- 
sion and lines from the primary not blended with the secondary. 
Then we varied the temperature of the secondary in order to repro- 
duce all identified lines of the secondary. We find that the observed 
spectrum is well fit with effective temperatures of 10500 ± 500 K 
and 5500 + 500 K for the primary and the secondary respectively, 
and with a ratio of radii equal to 1.5 ± 0.2, where the uncertain- 
ties are at the 2cr confidence level. For both components we used 
a surface gravity (g) given (in cgs units) by logg = 4.0, consis- 
tent with the evolutionary stage of these stars. The strong emission 
contamination does not allow us to measure this characteristic di- 
rectly. Our spectroscopic determination of the temperature of the 
primary is consistent with previous determinations of its spectral 
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Figure 2. Portion of the spectrum observed on August 26th 2005 (black 
full line) superimposed with the synthetic spectrum of the binary (red 
dashed line) calculated with a solar metallicity (upper) and a metaUicity 
of [M/H] = 0.5 (lower). The positions and the ions of the lines are indi- 
cated with vertical lines. The labels P and S beside the ion identification 
indicate lines of the primary and secondary, respectively. 



type, found to be between B9 and Al i Herbig & Kameswara Rad 



1 19721 : iFinkenzeller & Mundtl[T984l : Iffillenbrand et al.lll992h . 

We also find that the portion of the spectrum not contaminated 
with emission requires a higher metallicity than the sun in order 
to adequately fit the lines of the primary. The best-fit metallicity 
is [M/H] = 0.5. In Fig. |2l the upper panel shows the observed 
spectrum superimposed with a synthetic spectrum calculated us- 
ing BINMAGl and a solar metallicity, while in the lower panel 
the synthetic spectrum has been calculated using a metallicity of 
[M/H] = 0.5. It is however less evident that the secondary requires 
a high metallicity in order to fit its lines. While a higher metallicity 
seems to provide a marginally better fit to the unblended lines of 
the secondary, the differences between the solar metallicity and the 
[M/H] = 0.5 spectra are not sufficiently small as to be insignificant 
in comparison with the uncertainties of the observed spectrum. 

We can wonder if the high metallicity found in the primary 
reflects inaccurately determined parameters of the binary system. 
Using the ratio of radii and the temperatures of both components, 
we find that the derived luminosity ratio is equal to: = Lp/Ls = 
30t_f-j, which implies that the luminosity of the system is dominated 
by the luminosity of the primary. Therefore increasing Rp/Rg can- 
not result in a significant increase in the depth of the synthetic lines 
of the primary to provide a better fit to the observed lines. Further- 
more, we cannot sufficiently strengthen the computed lines of Fe, 
Ti and Si of the primary by changing the effective temperature or 
surface gravity within reasonable ranges. We therefore confidently 
conclude that the metallicity of V380 Ori A is higher than in the 
sun. 

It is well known that some magnetic intermediate mass stars 
on the main sequence show chemical peculiarities characterised 
by overabundances of elements such as iron, silicon and titanium. 
These stars are called the magnetic chemically peculiar Ap/Bp 
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Figure 3. (a) Ha profiles of the 24 observed spectra of V380 Ori. The date is indicated on the left side of the profiles. Note the temporal variations of the 
ampHtude, as well as the appearance and disappearance of substructures in the red and blue wings of the profile, (b) Difference between the Ha profiles and 
the mean Ha profile, (c) Same as (a) but plotted on a different scale to better show the absorption features in the blue wing of the profile. 



stars. First, as the primary of V380 Ori is magnetic, it is reason- 
able to consider that it might be a chemically peculiar star, or that it 
may evolve to become one. Secondly, the spectrum of V380 Ori A 
is dominated with iron, titanium and silicon lines, and the mea- 
surement of the metallicity described above has been mainly per- 
formed using these lines. The high metallicity that we measure may 
therefore reflect only overabundances in iron, titanium and silicon. 
Thirdly, it seems that the diffusion processes occurring at the sur- 
faces of chemically peculiar stars are sufficiently rapid to observe 
peculiar stars at the PMS stage of stellar evolution (Vink, Michaud 
and G. Alecian, private communication). However, at this stage the 
stars, still surrounded by gas and dust, might experience exchanges 
of matter with their surroundings, by accretion and winds, which 
could result in instabilities and in mixing of the stellar atmosphere. 
The diffusion processes are therefore unlikely to be efficient during 



the PMS phase. Indeed, we do not observe any chemical peculiari- 
ties at the surface of either of the youn g magnetic HAeBe stars HP 
200775 and HD 190073 jAlecian et al. 2008; C atala et al.ll2007h . 
whereas they are observed in the more evolved star HD 72106 (Fol- 
som et al. 2008) . The strong emission and variability observed in 
the spectrum of V380 Ori suggests that its environment is highly 
dynamic, with plenty of gas and dust. However, it is not clear if ex- 
changes of matter occur between the star and its environment. Fur- 
thermore the strong magnetic field that we observe at the stellar sur- 
face could sufficiently stabilise and isolate its atmosphere, to allow 
gravitational settling and radiative diffusion to occur Therefore, we 
tentatively conclude that the V380 Ori SB2 system is composed of 
a chemically peculiar primary star and a secondary of solar metal- 
licity. More observations of better S/N ratio, as well as a thorough 
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Figure 4. Portion of the spectrum observed on August 26th 2005 (black full 
line) and on February 12th 2006 (green dashed line). The positions and the 
ions of the lines are indicated with vertical lines. The label P beside the ion 
identification indicate lines arising from the primary. Note the variation in 
the emission of the metallic lines of the primary. 

study of the detailed chemical abundances of both components, are 
required to confirm this proposal. 

3.3 The emission lines 

3. 3. 1 Description 

The most evident characteristic of the Stokes / spectrum of V380 
Ori is a strong contamination by circumstellar emission. All Balmer 
and Paschen lines visible in the spectrum (from Ho- to He, and from 
PU to P20), as well as the Ca II 8498 A and 8662 A lines are filled 
with emission. Ho-, H/?, the Paschen lines and the IR Ca II lines ex- 
hibit a single-peaked emission profile. Ho- shows a faint blueshifted 
P Cygni absorption component that may originate from a wind (see 
Fig.[3](a) and (c) for Ho-). The radial velocity of the P Cygni com- 
ponent changes with time, from -350 kms"' to -250 kms"', with 
timescales of the order of few days. Our data reveal no periodicity 
of these variations. In Fig.[3](b) are plotted the differences between 
the Ho- profiles and the mean of these profiles, revealing changes 
in the shape and amplitude of the Ha profile on time scales of the 
order of days and sometimes hours (e.g. between both Ha profiles 
obtained on November 9th 2007 taken two hours apart). Complex 
structures seem to appear and disappear through the Ha profile on 
timescales of a day. The same kind of structures appear in the H/? 
profiles, in the same wing (either blue or red) as the Ha profile, but 
at different radial velocities. 

Most of the metallic photospheric lines in the spectrum of 
V380 Ori are also superimposed with emission. The emission is 
correlated with the intrinsic central depth and the excitation po- 
tential of the lines: at constant excitation potential, the greater the 
central depth, the greater the emission, while at constant central 
depth, the g r eater t he excitation potential, the lower the emission. 
ICatala et al.l ilOOt observe in the spectrum of another magnetic 
HAe star HD 190073, a similar correlation: the emission increases 
with the central depth of the lines. By re-analysing the spectrum 
of that star, we find also that at constant central depth: the greater 
the excitation potential, the lower the emission. These behaviours 
have also been observed in numerous non-magnetic HAeBe stars 
(Alecian et al. in prep.). The origin of these correlation is still not 
understood. 

Most Fe n and Ti ii lines in the spectrum show strong emis- 
sion with a presumably photospheric component in absorption, as 
illustrated in Fig. |4] The amplitude and the shape of these emis- 
sion lines changes on a timescale of days, and sometimes of hours. 




-400 -200 200 400 

velocity (km/s) 

Figure 5. Fe n 4351 A lines in the spectrum of V380 Ori observed on Au- 
gust 26th 2005 (black full line) superimposed with the fit (green dashed 
line). The 4 components of the fit, described in Sec. 3.3.1, are also plotted: 
the 3 Gaussian functions modelling the emission in blue dot-dot-dot-dashed 
lines, and the photospheric absorption in the red dot-dashed line. 



While the profiles of most of the lines are well reproduced with 
a simple Gaussian, the strongest emission lines require a slightly 
more complex model. These emission lines can be described as the 
superposition of three different Gaussian emissions: a strong main 
emission component superimposed with two smaller emission, one 
in each wing of the main component. We have fitted the strongest 
unblended emission lines using the following model. Each of the 3 
emission components is fit using a simple Gaussian with 3 free pa- 
rameters: the amplitude, the full-width at half-maximum (FWHM), 
and the centroid. When necessary, we also fitted the photospheric 
absorption using the convolution of a rotation profile (itself a func- 
tion o f V sin i and iVad) with a Gaussian of instrumental width l lGravl 
1 19921) . This convolution is then multiplied by a scale factor in or- 
der to fit the depth of the absorption component. An example of 
the results of this fitting procedure is showed in Fig. [5] While the 
FWHMs of the 2 wing Gaussians vary from one observation to an- 
other, and from one line to another, the strongest emission retains 
the same value of FWHM, around 65 kms"'. Temporal variations 
of the amplitude of the emission, and of the two wing emissions, 
are observed, but no periodicity is detected. The origin of these 
emissions is not known. The symmetry of the centroids of the wing 
emission with respect to radial velocity of the primary suggests that 
they may arise from a disk surrounding the central star. However 
the temporal variations of these emissions, i.e. variations of their 
amplitude, centroid, and FWHM, are difficult to explain in this con- 
text. While the main emission component of the iron and titanium 
lines is, most of the time, redshifted compared to the radial velocity 
of the primary star, we observe, in the same spectrum, a few lines 
for which the main emission component is centred or blueshifted 
compared to the radial velocity of the star. This behaviour is vari- 
able with time, and is not yet understood. Finally, in the few cases 
in which we include the photospheric component into the fit, the ra- 
dial velocities and the v sin / of this component are consistent with 
the values for the primary's photospheric lines found in Sec. 14.21 
by fitting the Stokes / LSD profiles. 

All the emission lines (except Ha, Hp, the Paschen lines and 
the IR Ca II lines) show a photospheric component, as already ob- 
served in the spectrum of HD 190073. Furthermore, the character- 
istic of the main emission component are similar to the emission 
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lines of HD 190073: in both cases the FWHM is around 65 km.s ', 
constant along the time and through the spectrum. As in Catala et 
al. (2007) we do not find a suitable scenario to explain these char- 
acteristics. These similarities cannot be due to a coincidence, and 
these emission properties bring therefore strong constraints on the 
forming region. 

Both Ca II H an d K lines sho w cl ear P Cygni profiles, 
as alre ady observed by Herbid d 19601) and iFinkenzeller & Mundtl 
( Il984) . The shape of the emission profile is also composed of a 
main emission component with wing components similar to those 
of the Fe ii and Ti ii lines. The radial velocity of the blue absorp- 
tion part of the P Cygni profiles also varies with time from -190 to 
-160 kms"' with respect to the rest wavelength of the Ca II lines. 

The He i 5876 A and He i 6678 A lines show broad emis- 
sion superimposed with a strong absorption component. Using a 
Gaussian function we have fitted only the emission component 
of these lines. We find that the equivalent width of the emission 
varies from 9 to 27 A for He i 5876 A, and from 10 to 18 A for 
He 1 6678 A, while their FWHMs vary from 130 to 200 km s"', and 
from 130 to 180 kms ', respectively. Their centroids vary from -30 
to 30 km s"', and are, most of the time, blueshifted with respect to 
the radial velocity of the primary. Interestingly, no emission, and 
only a small abso rption feature is visible at 587 6 A in the the spec- 
trum obtained by IFinkenzeller & Mundtl jl984l) . 

The O I 7773 A triplet is in emission, and superimposed with 
the photospheric absorptions. Its equivalent width varies with time 
from 3 to 7 A, with temporal variations similar to the He i and 
Balmer lines, but still without any periodicity. 
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Figure 6. HR diagram of V380 Ori. The red crosses are the error bars in 
temperature and luminosity of the primaiy (labelled P), and the secondary 
(labelled S). The PMS evolutionary tracks of 1, 1.5, 2, 3, and 4 Mq are plot- 
ted using full hnes. The zero-age main sequence (ZAMS; dot-dashed line), 
the isochrones (green dashed hnes) of 1 , 2 and 3 Myr, and the birthline (blue 
long dashed line), calc ulated with a mass accretion rate of 10"^ Mo.yr"' by 
IPalla & Stahled jl993l) . are also plotted. The filled diamonds indicate the 
moment where the mass of the convective envelope becomes lower than 
1 % of the stellar mass, and the triangles indicate the appearance of the con- 
vective core. 



3.3.2 Origin of the emission 

The angular separations between the primary and the two compan- 
ions (0.83 and 154 mas for the secondary and tertiary, respectively) 
are smaller than the pinhole aperture of the ESPaDOnS instrument 
(1.6")- Therefore all optical contributions from the 3 components 
are present in the spectra that we obtained. While no photospheric 
lines from the tertiary are detected in our spectra, we can won- 
der if a contribution from this component exists in the emission 
lines. iLeinert et al. IC993), using speckle interferometry, measured 
the flux of the primary and the tertiary separately in 6 IR photo- 
metric bands (V , J, H, K, L', and M). They fit the spectral energy 
distribution using two models: model a) is composed of a system in 
which the primary dominates the tertiary at all wavelength, whereas 
model b) assumes an optical primary and an IR companion. In both 
cases, the primary is dominant at optical wavelengths. It is therefore 
unlikely that the tertiary contributes to the emission in the optical 
spectrum of the system. 

Excluding the Balmer lines, the major emission in the spec- 
trum is observed only in the lines of the primary - such emission is 
a common characteristic of HAeBe stars. Our modelling of a few 
emission lines show that the wing emission components are sym- 
metric with respect to the radial velocity of the primary. Although 
the main emission component is not perfectly centred on the radial 
velocity of the primary, this characteristic is commonly observed in 
other HAeBe stars. Furthermore, the lines of the secondary that are 
not blended with lines from the primary show emission that is much 
weaker than that observed in the primary's lines. Finally, the vari- 
ability observed in the strong emission lines described above has 
been compared to the radial velocity variations of the secondary, 
and no correlation has been observed. The emissions observed in 
the metallic lines of the primary are therefore likely originating 
from the primary star of the system. However, the small separation 



(a sin ; < 0.33 AU) between both spectroscopic component, and the 
fact that the typical length of the disks surrounding HAeBe stars is 
larger than few AU ( Hillenbrand et al. 1992), lead to the conclusion 
that the disk would be circumbinary, instead of circumstellar. 

We cannot firmly conclude that this is true for the Balmer 
lines, and especially Ha and Hj8. The structure and the variability 
observed in these lines is complex and our observations cannot rule 
out that the secondary might contribute to them, an d even the ter - 
tiary as suggested by the spectro-astrometry study o f lBailevlll998h . 

The temporal variations observed among the emission features 
of the spectrum of V380 Ori are certainly the result of strong and 
rapid dynamic phenomena in the circumstellar environment of the 
star, which would be the subject of a very interesting study, but 
which is not the purpose of this paper. We therefore do not go fur- 
ther in the description and analysis of the emission features ob- 
served in the spectrum. Instead, we will look more closely at these 
characteristics in a future paper that will involve not only high S/N 
ratio spectroscopic observations, but also other type of observations 
such as linear polarimetry and interferometry, as well as sophisti- 
cated models capable of reproducing the spectroscopic emission of 
this star, and determining the physical parameters of its environ- 
ment. 

4 FUNDAMENTAL PARAMETERS OF THE PRIMARY 
AND THE SECONDARY 

4.1 The HR diagram 

V380 Ori belongs to the Ori OB lc association situated at a helio- 
centri c distance of 400 ± 40 pc faerbig & Bel]||l988l : iBrown et al] 
Il994h . The published photom etry shows that it is a variable star 
(V varying from 10.2 to 10.7. Ide Winter et alj[200lh . Its compan- 
ions and its environment are associated with significant dust and 
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Table 2. Fundamental parameters of the primary and secondary components 
of V380 Ori. Columns 2 to 5 give the effective temperature, the luminosity, 
the mass and the radius of the stars. Column 6 gives the age of the system. 
The projected rotational velocities (v sin i - derived in Sect. 4.2) are given in 
column 7. The macroturbulent velocity (Vma - also derived in Sect. 4.2) of 
the primary is indicated in column 8. v„,a of the secondary cannot be derived 
using our data. The 2tT uncertainties are indicated. 



Star Teff log(L/Lo) M/Mq RIRq age vsin; v^a 

(K) (Myr) (km.s"') (km.s"') 

P 10500 ±500 1.99 ±0.22 2.87^55^2 3 Q+i.i 6.7±1.1 3.5±0.8 

' 2±1 

S 5500 ±500 0.50^0 54 I.e+' O 2.0+'-^ 18.7 ±3.2 



gas, and this could be the cause of the observed variability. In this 
picture, the variability is due to variable extinction in the circum- 
stellar environment. Therefore, to compute the true luminosity of 
the V380 Ori system, we have adopted the bri ghtest magnitude 
in V and B, published bv Ide Winter et aL as the apparent 

magnitude of the system. This magnitude, in combination with a 
standard interstellar extinction {R\/ = 3.1), as well as the solar ab- 
solute magnitude My = 4.75, allow us to derive the luminosity 
of the system. We find log(L/Lo) - 2.00 + 0.20. From the lumi- 
nosity ratio determined in Sec. 3.2, we calculate the luminosity 
of the primary and the secondary: log(Lp/Lo) = 1.99 ± 0.22 and 
log(Ls/Lo) = 0.50!°;^^, respectively. 

Using these luminosities and the temperatures derived in 
Sec. 3.2, we compared the positions of both stars in the HR diagram 
(Fig.|6l with PMS stellar evolutionary tracks of sol ar metallicity, 
computed with the CESAM 2K code (' Morell [19971) , for dilferent 
stellar masses. From this comparison we infer the ages, masses and 
radii of both stars. The error bars in mass and radius are determined 
by taking ellipses around the error in temperature and luminosity. 
The age of the system is determined by taking the intersection of 
the age range of the primary (agep = 1.8^]^) and the secondary 
(agej = 7+^*), assuming that both stars have the same age. The 
results are summarised in Table|2] 

4.2 Calculation and fitting of the LSD Stokes / profiles 

In order to increase the signal to noise ratio of our data, to mea- 
sure accurate rotational and radial velocities of both stars, and to 
constrain their magnetic field character istics, we applied th e Least 
Square Deconvolution (LSD) method dOonati et alj|l997l) to our 
spectra using two different line masks (describing the positions, 
relative intensities and Lande factors of all lines predicted to be 
present in the spectrum), one for the primary component and one 
for the secondary. Both masks have been calcu lated using the stellar 
atmosphere code ATLAS9 of iKuruc j jl993h for the temperatures 
and gravities of both stars determined in Sect. 3.1. The masks have 
been cleaned as follows. For each star, we have selected lines that 
are not blended with the other component at any time, and that are 
not contaminated by emission, and we have rejected all other lines, 
including the Balmer lines. The resulting masks contain 85 and 374 
lines, respectively for the primary and the secondary component. 
The number of lines in the mask of the secondary is larger than in 
the primary mask due to its lower temperature, and also due to the 
fact that the spectrum of the secondary shows less contamination 
by emission. 

The results of this procedure were two sets of 24 Stokes /, 
V and diagnostic N LSD profiles, for the primary and secondary 



respectively. The peak signal to noise ratio in the Stokes V profiles, 
indicated in Table[T] are higher in the secondary profiles than in the 
primary's due to a larger number of lines in the secondary mask. 

To determine the radial velocity (Vrad) and the projected rota- 
tional velocity (v sin of both stars, we performed a simultaneous 
least-squares fit to the 24 LSD Stokes / profiles, independently for 
each star. Each profile was fit with the convolution of a rotation 
function (for which the projected rotational velocity v sin / is a free 
parameter in the fit) and a Gaussian whose width is c omputed fro m 
the spectral resolution and a macroturbulent velocity ( lGravlll992h . 

The free parameters of the fitting procedure are the centroids, 
depths, V sin i and macroturbulent velocity v^.^ of each component. 
The centroids of the functions can vary from one profile to an- 
other, whereas the depths, v sin i and Vn,a of each component cannot. 
This fitting procedure therefore assumes that the depths, v sin / and 
Vnia of both components do not vary with time, which we confirm 
(within the error bars) by fitting each profile separately. Including 
the macroturbulent velocity as a free parameter in the fit of the sec- 
ondary profiles results in a macroturbulence consistent with zero, 
but with a very large error bar. As the S/Ns of the secondary pro- 
files are very low, we therefore conclude that this parameter cannot 
be constrained by fitting these profiles. Therefore, for the fit of the 
secondary profiles, v^a has not been considered as a free parameter 
and has been fixed at km.s"' . 

Fig. |7] and [8] show the resulting fits superimposed on the ob- 
served Stokes / profiles of the primary and the secondary compo- 
nents, respectively. The profiles of November 6th and llth 2007 
have not been taken into account in the secondary fit because these 
profiles contain, in the red wing, a strong absorption component 
whose the origin is unknown. Including them in the fit results in a 
poorer fit and in larger error bars. We also observe that the depths 
of some profiles of the secondary are not well reproduced, which 
might be due to variable circumstellar emission or absorption con- 
taminating some spectral lines at few observational dates. However, 
the strongest differences observed between the fit and the observed 
profiles are within the error bars. Therefore, considering the com- 
plexity of the spectrum of the system, the mask used to compute 
these LSD profiles gives a satisfactory result. 

This automatic fitting procedure enables us to measure the 
vsin/ and radial velocities of both components. We obtain pro- 
jected rotational velocities of 6.7 ±1.1 and 18.7 + 3.2 km.s"' for the 
primary and the secondary components respectively, and a macro- 
turbulent velocity for the primary of 3.5 ± 0.8 km.s^l. The radial 
velocities of both components are included in Table[T] and are anal- 
ysed in the next section. 



4.3 Tlie radial velocity variations 

The separation between the primary and the tertiary being 200 
times larger than the separation between the primary and the sec- 
ondary, in the following, we will treat the system as a binary sys- 
tem, and we will ignore the impact from the IR companion in the 
radial velocity analysis. 

As mentioned in Sect. 3.3.1, the radial velocities of the spec- 
tral lines of the secondary vary with time, while those of the pri- 
mary do not. The radial velocities of the components, measured 
using the LSD Stokes / profiles, confirm this observation. We have 
fitted the temporal variation of the radial velocities of both compo- 
nents assuming a circular orbit. We should note that we also suc- 
ceeded in fitting these variations using an eccentric orbit. However, 
the uncertainties are sufficiently large that we have no indication 
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Figure 7. LSD Stokes / profiles of the primary stai' V380 Ori A, calculated from the 24 observed spectra (black full lines), superimposed with the best fit (red 
dashed line) described in Sec. 4.2. The dates of the observations are indicated beside each profile. 



I ' ' I 



**'|^^^^^^^^^^^j'^0^/2005 




26/8/2005 



9/1/2006 



/2006 



1 1/2/2006 



12/2/2006 



I ■ . I. ■ I ■ . ■ ■ I ■ . ■ ■ I 



l' ' I 




3/2/2006 



17/3/2007 



6/1 1/2007 



6/1 1/2007 



7/1 1/2007 



7/11/2007 



. I. . I 



l' ' I 




8/1 1/2007 



8/1 1/2007 



9/1 1/2007 



9/11/2007 



10/1 1/2007 



10/11/2007 



I . . . . I . . I. . I 




1/1 1/2007 



1 1/1 1/2007 



2/1 1/2007 



12/1 1/2007 



13/1 1/2007 



13/1 1/2007 



-50 50 100 150-50 50 100 150-50 50 100 150-50 50 100 150 



veiocilj' (i<in/s) 



velocity (km/s) 



velocity (i<in/s} 



velocity (km/s) 



Figure 8. LSD Stokes / profiles of the secondary star V380 Oii B, calculated from the 24 observed spectra (black full fines), superimposed with the best fit 
(red dashed line) described in Sec. 4.2. The dates of the observations are indicated beside each profile. The profiles of November 6th and November 1 1th 2007 
show an absorption component in the red wing, whose the origin is unknown. Therefore these profiles have not been included in the fit. 



that these orbits are eccentric, and we will only consider the sim- 
plest model that can reproduce our measurements. 

As the radial velocity of the primary does not vary signifi- 
cantly, we chose to include the measurements into the fit as a con- 
stant with time equal to the systemic radial velocity, in order to 
better constrain this parameter. We performed a simultaneous fit of 
both radial velocity curves using simple sinusoids for each curve. 
This fitting procedure depends on 4 parameters: the orbital period 
(Porb) of the system, the reference time (To), the systemic radial ve- 
locity (y), and the radial velocity amplitude of the secondary (A's). 
We determined an upper limit on the radial velocity amplitude of 



the primary (Kp) using the mean of the error bars of the measure- 
ment of its radial velocity. 

The radial velocities of both components, and the result of 
the fit of the orbits, are plotted in Fig. |9l folded with the orbital 
period. The value of the fitting parameters are presented in Table 
[3] We derive only a lower limit on the mass ratio Mp/M^ of 2.9, 
consistent, within the error bars, with the mass ratio obtained from 
the HR diagram in Sec. 4.1 (Mp/Ms = l-S^^i g)- Using the masses 
derived from the HR diagram, and the inferred value of Mp sin^ 
and Ms sin' /c-i,, we find that the inclination angle between the orbit 
axis and the observer's line of site (i'orb) is lower that 31°, using the 
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Figure 9. Radial velocities of the primary (blue points) and the secondary 
(red triangle) plotted as a function of the orbital phase (according to the 
ephemeris given in Table[3), and superimposed with the best-fit radial ve- 
locity curve of the secondary. The systemic radial velocity is indicated with 
a black dashed lines. The lower panel shows the radial velocities of the 
primary on a smaller scale. 



Table 3. Orbital parameters of the system, a is the distance between both 
stars, ap and as are the radii of the primary and the secondary (circular) 
orbits, respectively, /ori, is the inclination angle between the orbit axis and 
the line of site of the observer. has not been determined within the fitting 
procedure, but from the error bars of the radial velocities of the primary. The 
3iT error bars are indicated. 





Results of the fit 


Period (days) 


104 ±5 


To (HDJ) 


2453505 ± 18 


r (km/s) 


27.8 ±2.1 


Ks. (km/s) 


18 ± 14 


Kf (km/s) 


< 1.4 




Derived parameters 


Mp/Ms 


>2.9 


a sin i'opi, (AU) 


< 0.33 


ap sin lopb (AU) 


< 0.0139 


as sin ('orb (AU) 


0.17 ±0.14 


Mp sin^ !ojb (A/o) 


< 0.40 


Ms sin^ i„h {Mq) 


< 0.018 



primary mass, and lower than 13°, using the secondary mass, which 
indicates that the orbit might be seen close to pole-on. Due to the 
relatively low signal to noise ratio of our data, the uncertainties 
on the orbits and fundamental parameters of the stars are large. 
More spectroscopic observations of higher signal to noise ratio are 
required in order to better constrain the orbital parameters and the 
mass ratio of the system. 
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Figure 10. LSD Stokes /, V, and diagnostic N profiles of the primary (left) 
and the secondary (right), normalised to the continuum of the binary 7^. 
These profiles have been obtained from the first November 12th 2007 spec- 
tram. The mean error bars are indicated beside each profile. The V and N 
profiles have been shifted vertically and have been multiplied by a factor of 
3 for the primary and 2 for the secondary, for display purpose. 



5 MAGNETIC FIELD ANALYSIS 

Fig. [To] shows the LSD Stokes /, V, and diagnostic profiles of 
the primary (left) and the secondary (right) for one night, calcu- 
lated using the masks described in Sec. 14.21 We observe that the 
null A' profile is totally flat in both cases, indicating that the signa- 
ture observed in the Stokes V profile of the primary is real, and is 
not a spurious signal created by instrumental artefacts. The mag- 
netic field detection in the primary is furthermore confirmed by the 
detection of Zeeman si gnatures in the H o- line and in each line of 
the triplet 01 777 nm ( Wade et al. 2005). Unlike the primary, the 
Stokes V profile of the secondary doesn't show any Zeeman signa- 
ture, meaning that no magnetic field has been detected at the surface 
of the star. It also confirms that the mask of the secondary does not 
contain lines contaminated with lines from the primary. 

In order to characterise the magnetic field of the primary and 
to put an upper limit on the strength of the magnetic field that the 
secondary could host, we must perform an analysis of the Stokes 
profiles of both stars, normalised to their individual continua. The 
observed spectrum is the combined light coming from both stars, 
which means that the intensity in the continuum (7^) is the sum 
of the continuum intensity of both stars: 4 = 4? + hs,. The lu- 
minosity of the system being dominated by the luminosity of the 
primary, the contribution from the secondary's continuum is neg- 
ligible compared to the error bars in the Stokes / and V profiles 
of the primary, and won't be taken into account in the analysis of 
these profiles. However, we cannot neglect the contribution from 
the primary in the Stokes profiles of the secondary. We have there- 
fore re-normalised these profiles, to the continuum of the secondary 
only, using the luminosity ratio determined in Sec. 13.21 As we have 
been careful to not include lines blended with the companion in the 
calculation of the LSD profiles, the intensity inside the lines of the 
secondary is only contaminated with the continuum (and not the 
lines) of the primary. Therefore the intensity and circular polarisa- 
tion inside the lines of the secondary are given by: 

!_ ^ 4p + Is ^ r^ + hlhs . V Vs Vs 1 



and — 



hp + /cs 1 + a h hp + hs hs 1 + rL 

where hp and hs are the continuum intensity of the primary and 
the secondary respectively, /s is the intensity from the secondary 
only, Ks is the Stokes V profile of the secondary, and is the lumi- 
nosity ratio as defined in Sec. 13.21 Using these equations, the LSD 
profiles calculated from the combined spectrum, and the value of 
previously determined, we can extract the Stokes profiles of the 
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Figure 11. Longitudinal magnetic field {Be) of the primary plotted as a 
function of rotation phase, according to the ephemeris summarised in Sect. 
5.1. The best fit sinusoid is also plotted as a red dashed line. The data points 
for both November 9th 2007 observations (with relatively large eiTor bars) 
have been removed from the plot, for display purpose. 

secondary normalised to the continuum of the secondary: /s//cs and 
Vs//cS- In the following we describe the analysis performed on the 
profiles of both stars. 



5.1 Phased longitudinal field variation of the primary 

In order to characterise the magnetic field of the primary, we need 
to model the Stokes / and V profiles. The magnetic fields of the 
Ap/Bp stars and of the other field HAeBe stars are large-scale or- 
ganised fields, with an important dipole component. We therefore 
used the oblique rotator model, i.e. a dipole of polar strength Sd 
placed at the centre of the star, with a magnetic axis inclined at an 
angle fi with respect to the rotation axis, and the line of sight mak- 
ing an angle i with the rotation axis (Landstreet 1970). As the star 
rotates with a period P, the magnetic configuration on the observed 
surface of the star changes with time. Therefore the Stokes V pro- 
file, and the magnetic strength projected onto the line of sight and 
integrated over the observed surface, i.e. the longitudinal magnetic 
field (Sf), change with time. 

As described in Sect. 2, we monitored V380 Ori during many 
nights spread over about 3 years, in order to fully sample the ro- 
tational cycle of the star. We measure d the longitudinal magnetic 
field of the 24 profiles using Eq. (1) of lWade et af] j200(]t) : 

/ V{v)vdv 

Sf = -2.14 X 10'^ ^ , (1) 

g^c/(l-/(v)//,)dv 

where the wavelength A in A, and the Lande factor g are the mean 
values calculated using all the lines of the mask used to compute 
the LSD profiles, Be is in gauss, and c is the speed of light. We 
calculated the first moment of the Stokes V profile and the equiva- 
lent width of the intensity profile, with an integration range from 12 
to 45 km s"' . We found values of Bi ranging from -600 to 700 G 
with uncertainties between 60 and 300 G (see Table[TJ. Using &x^' 
minimisation technique, we fit the variation of Be with a sinusoidal 
function of four parameters: 

• P, the rotation period of the star, 

• fo, the reference Julian Date of maximum of Be, 

• 5, the semi-amplitude of the curve. 



• Bo, the vertical shift of the sinusoidal curve with respect to 
Sf = 0, 

consistent with the oblique rotator model. The best-fit model, su- 
perimposed on the data in Fig.[TT] gives the following parameters: 
P = 4.43 ± 0.23 d, to = 2453413.40 ± 0.30 d, B = 440 ± 200 G, 
So = 0± 140 G, with^^ = 1.4. In contrast, the^^ calculated assum- 
ing a zero-field model (in which the longitudinal magnetic field is 
constant and equal to 0) is 1 1.3. 

The detailed measurements of the longitudinal magnetic field 
values depend sensitively on the integration range used in the cal- 
culation of the first moment of V and the equivalent width. We have 
measured Be for numerous integration ranges, by varying the blue 
limit from to 18 km,s ', and the red limit from 40 to 60 km.s"'. 
We find that the standard deviation of the values of Be measured 
from a single LSD profile can be as high as 285 G. Furthermore, 
the measured values for two profiles observed on the same night 
can show important differences as large as 700 G, while the Stokes 

V profiles do not show any significant difference. These effects can 
have important consequences for the determination of the ampli- 
tude of variation B, as well as the shift of the sinusoid Bo- In con- 
trast, however, whatever the integration range used, the Be values 
show similar variations with a periodicity around 4.4 d, implying 
that the determination of the period is insensitive to the integra- 
tion range. These characteristics have already been observed dur- 
ing our previous an alysis of the magnetic HBe star HD 200775 
jAlecian et al.l2008h . As for HD 200775, we won't try to determine 
the magnetic geometry of V380 Ori A using the longitudinal mag- 
netic field, and we will determine it by fitting directly the Stokes 

V variations, using the period derived in this fit as a first guess in 
the Stokes V modelling. We should note that Fig. [TT] shows that 
the phase coverage is uniform and complete, which is optimal for 
modelling the magnetic topology of the star. 



5.2 Fitting of the Stokes V profiles of the primary 

The observed Stokes V profiles are plotted in black full lines in Fig. 
[T2I We see clear variations indicating rotational modulation of the 
magnetic field at the visible surface of the star. To fit the Stokes V 
profiles directly, we use the dipole oblique rot ator model described 
above. Using the relations of lLandstreej jl97ol) . giving the intensity 
of the magnetic field at each point on the stellar surface, we calcu- 
late the longitudinal magnetic field bi{6, <p), at each point (6,(fi) of 
the surface of the star, in classical spherical coordinates in the ob- 
server's frame. 

We assume a gaussian local intensity profile of width cr 
and depth d. The width is calculated using the resolving power 
of the instrument and the macroturbulent velocity determined 
in Sec. 14.21 The depth is determined by fitting the Stokes / 
profiles of the primary component, determined in Sec. 14.21 We 
then calculated the local Stokes V profile at each point on the 
surfa ce of the star using the weak magnetic field relation of 
Landi degl'Innocenti & Landi de gl'In nocenti, Cl973i) : 

dl(e, w) 

V(v, 0, v) = -CgAocbe{9, f)^T^, (2) 
dv 

where C = 4.67 x lO^'^A^'G"', ^ and Aq (A) are the mean Lande 
factor and wavelength of the lines used in the mask (Section 2), 
c is the speed of light, 1(8, ip) is the local intensity profile, and v 
is the local radial velocity due to stellar rotation. Then we inte- 
grated over the visible stellar surface usin g the limb darkening law 
with a parameter equal to 0.5 ( IClaretl2000l) . We obtain the synthetic 
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Figure 12. LSD V profiles of the primary (full black lines) superimposed by the best oblique rotator model (dashed red line). The numbers close to the profiles 
ai'e the rotation phase, and the small bars on the right of the profiles are the mean error bars in V. The profiles ai'e sorted by increasing rotational phase, and 
the date of observation is indicated next to each profile. 



Stokes V(i') profile, that we normalise to the intensity continuum, 
to compare to the observed V profiles. This model depends on five 
parameters : 



Table 4. Ranges and minimum bins of parameters explored in the fit of the 
Stokes V profiles 



• P the rotation period 


Parameters 


Min 


Max 


bin 


• to, the reference Julian Date of the maximum of the surface 


P (days) 


3.0 


5.0 


0.00001 


magnetic intensity, used with P to compute the rotation phase 


To (HJD) 


2453510 


2453515 


0.001 


• (, the inclination of the stellar rotation axis to the observers 


in 





180 


0.1 


line-of-sight. 







180 


0.1 


• y6 the magnetic obliquity angle. 


Bd(G) 





6000 


10 



the dipole polar magnetic intensity. 

We calculated a grid of V profiles for each date of observation 
(see Table |4] for details on the grid), varying the five parameters 
and assuming for the initial value of P the solution obtained from 
modelling the longitudinal field variation. Then we applied a x' 
minimisation to find the model which best matches simultaneously 
our 24 observed profiles. The best model that we found, ^i\hx~ = 
1.6, corresponds to P = 4.31276 + 0.00042 d, i = 32 ± 5°, /3 = 
66 ± 5, and = 2120 ± 150 G, where the error bars correspond 
to the 3cr confidence level. Figure[T2]shows the synthetic Stokes V 
profiles superimposed on the observed ones. We see that this model 
acceptably reproduces most of the observed V profiles. 

The period found using this fitting procedure is consistent with 
that derived from the longitudinal magnetic field variation. How- 
ever, neither / nor /J is close to 90°, a result which is inconsistent 
with the symmetry of the longitudinal magnetic field variation with 
respect to Bf = . However, the amplitude of the Stokes V pro- 
files around phases and 0.5 are clearly different, leading to differ- 
ent absolute values of the longitudinal magnetic fields. From these 
simple observations we therefore do not expect a symmetric lon- 
gitudinal magnetic field curve. We are therefore confident with the 
magnetic configuration found using the fitting of the Stokes V pro- 
files, and we confirm that using Bi values obtained from low S/N 
ratio data are unreliable for determining the detailed magnetic con- 
figuration of a star. 



5.3 Stokes V analysis of the secondary 

We used the Stokes / and V profiles of the secondary star to de- 
termine the upper limit on the strength of a magnetic field that the 



Table 5. Magnetic dipole model of V380 Ori A. The 3cr error bars are 
indicated. 



P (days) 
To (HJD) 

in 



4.31276 ±0.00042 
2454412.997 ±0.043 
32 ±5 
66 ±5 
2120 ± 150 



star could host. According to the position of the secondary in the 
HR diagram (Fig. |6l(, its PMS evolutionary status is not well de- 
fined. The star might possess a convective envelope under the sur- 
face that could give rise to magnetic fields of complex structure, 
through a convective dynamo as in the T Tauri stars. Therefore we 
do not claim that this star hosts a large-scale organized magnetic 
field, however as no magnetic field has been detected we will as- 
sume the simplest model, which is the oblique rotator model as 
described above. 

Our method consists of calculating simulated Stokes V obser- 
vations of the same S/N ratio as our data. For each one of the 24 
LSD profile sets, and for a fixed value of the dipole strength Bj, 
we have calculated 100 simulated Stokes V profiles as follows. For 
each of the 100 iterations, we adopted random values of ; and [}. We 
then computed 24 synthetic Stokes V profiles using the oblique ro- 
tator model, assigning a random phase to each LSD profile. Jn this 
calculation we used a Stokes / profile computed with the vsin/, 
central depth and macrotubulent velocity found in Sec. 14.21 a radial 
velocity equa l to km s"' , and a limb-darkening coefficient of 0.7 
( IClarelll200(ll) . To the computed Stokes V profiles we then added 
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synthetic Gaussian noise consistent witli tlie S/N ratio of the ob- 
servations. For each of these simulated profiles we calculated the 
detection probability (Donati et al. 1997). We consider that a mag- 
netic field would have been detected at 3(T if the detection probabil- 
ity is greater than 0.9973. We calculated the number of iterations 
that gives at least one Scr detection among the 24 simulated pro- 
files, and we repeated this procedure for difi"erent values of B^. We 
find that if is greater than 800 G, the probability to have at least 
one detection among our data is 99%, while if B^ is greater than 
500 G this probability is 90%. Therefore, if the secondary host a 
dipole at its centre, its strength is very likely to be below 500 G. 



6 DISCUSSION AND CONCLUSIONS 

This paper describes the analysis of high-resolution spectropolari- 
metric observations of the pre-main sequence Herbig Ae star V380 
Ori, obtained with the new generation instruments ESPaDOnS and 
Narval. The main goal of this stud y was to char acterise the stellar 
magnetic field discovered in 2005 tWade et al.l2005i) . and to further 
explore the origin of the magnetic fields of intermediate mass stars. 
The results of the analysis are summarised and discussed below. 

First, a careful investigation of the intensity spectrum of V380 
Ori led to the discovery of a cool spectroscopic companion, and 
confirmed the detection of the Li i 6707 A line reported by previ- 
ous authors. We argue that this spectro scopic companion ca nnot be 
the infra-red companion discoverd bv lLeinert et al. which 
implies that V380 Ori is a triple system with two spectroscopic 
components and an IR companion. 

The modelling of the binary spectrum allowed the determi- 
nation of the effective temperatures of both stars, as well as the 
luminosity ratio of the system, and revealed a high metallicity 
([M//f] = 0.5) of the primary. The same metallicity was also able 
to reproduce the lines of the secondary, but was not required to ob- 
tain an acceptable fit of its spectrum. As a magnetic field has been 
discovered only in the primary, we therefore proposed that the pri- 
mary is a magnetic chemically peculiar star, as observed among 
the magnetic A/B star on the main sequence. No variability is ob- 
served in the shape and strength of the lines, implying that any 
non-uniformities in the surface distribution of Fe and Ti are ei- 
ther too small to be detected, or that these elements are distributed 
symmetrically about the rotation axis. These results, if confirmed 
with more accurate observations and a thorough study of the abun- 
dances, would give the first Herbig Ae star, whose the environment 
is highly dense, that is able to develop abundance anomalies at its 
surface. It would also lead to unique constraints on the physical 
processes at the origin of these anomalies and the timescales for 
developing them, as well as constraints on the interaction of the 
star with its environment in the presence of a strong magnetic field. 

The intensity spectrum was found to be strongly contaminated 
with emission. The emission in the metallic lines is correlated with 
the central depth and the excitation potential of the lines: at constant 
excitation potential, the greater the central depth, the greater the 
emission, while at constant central depth, the greater the excitation 
potential, the lower the emission. The strongest emission, observed 
mainly in the Fe n and Ti ii lines, shows 3 sub-structures, each of 
which was well reproduced with a Gaussian function. The fit of 
these lines revealed 

• a main emission component with a constant FWHM of 
~65 km.s"' over the spectrum and over time, slightly redshifted 
compared to the radial velocity of the primary. 



• two wing emissions, one in each wing of the main emission 
component, whose amplitude and broadening vary with time, that 
might come from a disk surrounding the central star. 

We argued that these emissions are very likely produced in a 
circumbinary disk whose the emission properties are mainly driven 
by the primary's light, and might not be contaminated with the 
emission of the secondary. However, the three components of the 
system might contribute to the strong emission profiles of com- 
plex structure and variability observed in the Balmer lines. Our data 
did not allow us to determine any correlation between the temporal 
variations of the wing emission components and the rotation period 
of the primary star. We are therefore unable to assess whether the 
proposed disk surrounding th e primary is produced b y magnetically 
confinement, as proposed bv lTownsend et al] ( l2005h to explain the 
disk surrounding the hotter main sequence Bp star a Ori E. On the 
other hand, the lack of veiling (associated with accretion shock hot 
spots), the lack of correlation between the rotation period and emis- 
sion variability, and the absence of distortion in the Stokes / profiles 
(signs of hot or cold spots) suggest that magnetospheric accretion is 
not occurring around the primary star, as one might expect from its 
PMS status and the similarity between its emission spectrum and 
those of the T Tauri stars. 

The strong similarity of the characteristics of the emission 
lines between V380 Ori and HD 190073 is puzzling. If the main 
emission component is arising in a disk, we would expect to mea- 
sure different FWHM of their emission profile, due to different in- 
clination angle. The probability to observe two magnetic HAe stars 
under the same angle being very low, the main emission compo- 
nent might come from another forming region with isotropic prop- 
erties. As in Catala et al. (2007) we can consider the scenario of a 
heated region at the base of a wind arising from the photosphere, 
with turbulent motions. However, the problem encountered with 
HD 190073 is also encountered with V380 Ori: most of the emis- 
sion lines are redshifted with respect to the radial velocity of the 
star, while we would expect them to be blueshifted if they come 
from a wind. The similarity of the emission spectrum of both stars 
therefore puts very strong constraints on the physical and geomet- 
rical properties of their circumstellar environment. 

The P Cygni absorption observed in Ho- and in the Ca H and K 
lines revealed the presence of a wind surrounding either the central 
star, or the central binary. The shape, radial velocity and broaden- 
ing of all these emissions showed strong variations on time-scales 
of days, and sometimes hours, suggesting a vigorous dynamic in 
the environment of the stars. More observations of better S/N ra- 
tio coupled with linear polarimetric and interferometric data, and 
compared to models of circumstellar disks and winds, would be re- 
quired in order to conclude about the circumstellar geometry and 
the physical processes involved in the formation and the dynamic 
of these emissions. 

Using the LSD Stokes / profiles of both components we de- 
termined their projected rotational velocities and radial velocities. 
We modelled the variations of the radial velocities of both stars and 
determined the orbital parameters and the mass ratio of the system. 
As no variations are observed in the radial velocity of the primary, 
only a lower limit was determined on the mass ratio, a value that is 
consistent with the mass ratio determined using the HR diagram. 

The LSD Stokes V profiles of the primary show a clear Zee- 
man signature characteristic of an organized magnetic field at the 
surface of the star. We found that the oblique rotator model of a 
centred dipole was able to reproduce acceptably the temporal vari- 
ations of the Stokes V profiles, implying that the large-scale mag- 
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netic field of V380 Ori A is a simple dipole, as observed in the other 
Herbig Ae/Be stars and in the Ap/Bp stars. The result reinforces ar- 
guments in favour of a fossil magnetic origin by showing that there 
is no exception to the apparent rule that all magnetic intermediate 
mass stars host large-scale, organised magnetic fields. 

The fit of the Stokes V profiles led to an inclination of the 
rotation axis equal to 32 ± 5°. However, assuming rigid rotation 
and using the derived rotation period of 4.31276 ± 0.00042 d, the 
radius inferred from the HR diagram position, and the v sin i mea- 
sured from the LSD Stokes / profiles, we infer i = 1 1+^°, which is 
inconsistent with the determination from the Stokes V fit. In order 
to understand this difference, we performed difi"erent tests. First, 
we fixed the inclination angle to 11°, and we performed a new fit 
on the Stokes V profiles. The result produced synthetic Stokes V 
profiles inconsistent with the observed ones around phase 0.0. Sec- 
ondly, the plot of the as a function of the period shows several 
local minima, and therefore we searched for the best magnetic con- 
figuration, reproducing the Stokes V profiles around the periods as- 
sociated with the local minima. In each cases, the synthetic profiles 
were found to be inconsistent with the observed profiles. However, 
the magnetic strength inferred when we performed these tests re- 
mained very close to our final determination (Sj = 2100 G). We 
therefore conclude that the magnetic configuration and the rotation 
period that we find using the fit of the Stokes V profiles is robust 
and is not the source of the disagreement between both values of 
the inclination angle. 

This disagreement might comes from the determination of the 
radius of the primary, inferred from the effective temperature and 
the luminosity of the star. While we are confident with the determi- 
nation of the temperature within the error bars, our determination 
of the luminosity is more debateable. The error bars on the lumi- 
nosity are determined from the uncertainties on the distance and 
the luminosity ratio (Lp/Ls). We did not take into account the fact 
that the star is variable, as well as the uncertainty on the interstel- 
lar extinction, which might lead to additional uncertainties. We en- 
countered the same problem during our study of HD 200775: an 
inco nsistency betw een two different determinations of the inclina- 
tion jAlecian et ^HoOS) . As for HD 200775, we must conclude 
that the determination of the luminosity of HAeBe stars, as well 
as the fundamental parameters using the HR diagram, are highly 
uncertain, and must be used with caution. 

If the luminosity that we determined is overestimated, the in- 
ferred mass and radius of the primary might be also overestimated. 
Using our estimation of the mass and using the PMS stellar evolu- 
tionary models, we can determine the radius of the primary when it 
will reach the ZAMS. Notwithstanding that our estimates of masses 
and radii of the components are highly uncertain, we can reason- 
ably assume that the ZAMS radius of a 10500 K star is greater than 
2 Rq. Using our determination of the radius, and the fact that we 
might have overestimated it, the upper limit on the actual radius is 
4 Rq. We therefore deduce that the ratio of the actual radius and the 
ZAMS radius is not larger than 2, and can be as low as 1 . Using 
this range of R/Rzams, we can derive the magnetic strength and the 
period on the ZAMS. 

Assuming magnetic flux conservation during the PMS phase, 
we predict that, when the primary will reach the main sequence, 
its magnetic strength will be larger than ~ 2 kG, which is of the 
same order of magnitude as the magnetic strengths of the Ap/Bp 
stars, as expected in the context of the fossil model. Assuming the 
conservation of angular momentum, we find that the rotation period 
on the ZAMS can be as low as 1 . 1 d, similar to the rotation periods 



of the Ap/Bp stars (> 1 d), meaning that V380 Ori A might evolve 
towards the ZAMS with a constant angular momentum. 

Assuming that / = 32°, the true equatorial velocity of the pri- 
mary is estimated to be 13 kms"', showing that V380 Ori A is 
already a slow rotator at an age around 2 Myr. In order to e xplain 
the slow rotation observed among the Ap/Bp stars, IStepienI Owdti 
proposed that magnetic braking occurs during the PMS phase of 
stellar evolution. As V380 Ori A appears already to have lost sig- 
nificant rotational angular momentum, our results show that if this 
star experienced a braking mechanism it must have been during the 
very first stages of the PMS evolution, or before. 

Using simulated observations of Stokes V profiles of the sec- 
ondary, calculated from the oblique rotator model, we determined 
an upper limit of 500 G (at 90% confidence) on any dipolar mag- 
netic field that this star may host. If the fossil model is correct, 
one would expect that a binary system should be composed of two 
magnetic stars with similar magnetic flux, if both stars have formed 
from the same molecular cloud (which is a reasonable hypothesis 
considering the young age of this system). Using the radii deter- 
mined in Sect. l4.1l we find an upper limit of 145 kG R^ on the mag- 
netic flux (B^xR^) of the secondary, and a lower limit of 100 kGR^ 
for the magnetic flux of the primary, which means that our results 
are consistent with the fossil field hypothesis. However, the radii 
of the stars are very poorly constrained. In order to better constrain 
the fossil model we need to improve the determination of the funda- 
mental parameters of both stars by performing interferometric ob- 
servations, and measuring separately both luminosities. However, 
the instrumentation available today is limited to bright magnitudes 
and won't be able to resolve the stars. In few years, once the effi- 
ciency of the interferometers has improved, these observations will 
be feasible and could be leaveraged to bring new constraints to test 
the fossil hypothesis. 
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